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ABSTRACT 

Far-infrared observations provide crucial data for the investigation and characterisation of the 
properties of dusty material in the Interstellar Medium (ISM) , since most of its energy is emitted 
between ~ 100 and 200 /Ltm. We present the first all-sky image from a sensitive all-sky survey 
using the Japanese AKARI satellite, in the wavelength range 50 - 180 /xm. Covering > 99% of 
the sky in four photometric bands with four filters centred at 65 /um, 90 fiva, 140 fim, and 160 /im 
wavelengths, this achieved spatial resolutions from 1 to 2 arcmin and a detection limit of < 10 
MJy sr _1 , with absolute and relative photometric accuracies of < 20%. 

All-sky images of the Galactic dust continuum emission enable astronomers to map the large-scale 
distribution of the diffuse ISM cirrus, to study its thermal dust temperature, emissivity and 
column density, and to measure the interaction of the Galactic radiation field and embedded 
objects with the surrounding ISM. In addition to the point source population of stars, protostars, 
star-forming regions, and galaxies, the high Galactic latitude sky is shown to be covered with 
a diffuse filamentary-web of dusty emission that traces the potential sites of high latitude star 
formation. We show that the temperature of dust particles in thermal equilibrium with the 
ambient interstellar radiation field can be estimated by using 90 fim, 140 fim, and 160 fim 
data. The FIR AKARI full-sky maps provide a rich new data set within which astronomers can 
investigate the distribution of interstellar matter throughout our Galaxy, and beyond. 
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1. INTRODUCTION 

Infrared continuum emission is ubiquitously observed 
across the celestial sky and is attributed to the thermal 
emission from interstellar dust particles. The spectral 



energy distribution (SED) of the dust continuum emis- 
sion has been observed for various astronomical objects 
including diffuse interstellar medium, star-forming re- 
gions and galaxies. The dust continuum SED has its 
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peak at far-infrared (30-300 /im; FIR) wavelengths, 
with approximately two thirds of the energy being ra- 
diated at A > 50 /im (Draine 2003; also see Compiegne 
et al. 2011). Since the stellar radiation is mainly ab- 
sorbed by large grains (LGs) in interstellar space, and 
converted to FIR continuum emission, it is important 
to measure the total FIR continuum emission energy 
as it is a good tracer of total stellar radiation energy 
and thus a good indicator of the star-formation activ- 
ity (Kennicutt 1998). The = 2 modified black-body 
spectrum of temperature 15 - 30 K has its peak at 100 
- 200 /im. It is thus essential to observe the FIR con- 
tinuum, including this wavelength range, to measure 
the total FIR emission energy from LGs to trace the 
star-formation activity. 

An all-sky survey of FIR continuum was pioneered 
by the IRAS satellite (Neugebauer et al. 1984), which 
was, until the completion of the AKARI mission, the 
only available FIR all-sky survey with high spatial res- 
olution. Given the importance of the FIR photometric 
data, it is valuable to upgrade the data by improv- 
ing spatial resolution, and by adding more photometric 
bands, especially at longer wavelengths. 

For that purpose, we performed new all-sky survey 
observations with A KARI (Murakami et al. 2007). In 
this paper, we describe the quality of the image data 
and its capability to investigate spatial distribution of 
ISM and its SED with improved spatial resolutions for 
the whole sky. 

2. Observation and Data Reduction 

The AKARI satellite (Murakami et al. 2007) was 
launched in February 2006 and the all-sky survey was 
performed between March 2006 - August 2007, which 
was the cold operational phase of the satellite with 
liquid Helium cryogen. The all-sky survey observa- 
tions in FIR wavebands were carried out by using an 
FIR instrument called the Far-Infrared Surveyor (FIS: 
Kawada et al. 2007), which covered the 50 - 180 //m 
waveband (Kawada et al. 2007). 

Four photometric bands are used for the all-sky sur- 
vey. Two of four bands had a broader wavelength cov- 
erage and continuously covered the whole waveband 
range: the WIDE-S band (50-110 /im, centred at 90 
/im) and the WIDE-L band (110-180 /im, centred at 
140 /xm). The other two bands had narrower wave- 
length coverage and measured both the shorter and the 
longer ends of the wavebands: the N60 band (50-80 



Table 1. 

Scan coverage of the survey observations 



Spatial scans Single 


Multiple 


> 5 


> 10 


Unscanned 


N60 


99.1% 


96.9% 


60.3% 


13.5% 


0.9% 


WlDE-S 


99.1% 


97.0% 


61.0% 


13.8% 


0.9% 


WiDE-L 


99.5% 


98.4% 


78.7% 


25.9% 


0.5% 


N160 


. . . 99.5% 


98.4% 


76.8% 


24.1% 


0.5% 


4 BANDS . . . 


. . . 99.1% 


96.8% 


60.0% 


13.2% 


0.9% 



/im, centred at 65 /im) and the N160 band (140-180 
/im, centred at 160 /im). 

The observed fraction of the sky during the survey 
period (the cold operational phase) is tabulated in Ta- 
ble 1. About 97 % of the sky was multiply surveyed in 
the four photometric bands. 

The data obtained during the observations was pre- 
processed using the AKARI pipeline tool, which was 
originally optimised for point source extraction (Yama- 
mura et al. 2009). Following this initial pre-processing, 
a separate pipeline was then used to recover the large- 
scale spatial structures accurately. The main proper- 
ties of the diffuse mapping pipeline are: 1) transient 
response correction of the detector signal, 2) subtrac- 
tion of zodiacal emission, 3) image processing, 4) image 
destriping, and 5) recursive deglitching and flat-fielding 
of processed image. Processes 1) and 2) were performed 
on the time-series data and 3) - 5) were performed on 
the image plane data. The details of these procedures 
are found in Doi et al. (2009) and Doi et al. (2012). 

3. Results 

Intensity maps from the all-sky survey in the four pho- 
tometric bands are shown in Fig. [TJ Detailed photom- 
etry with high spatial resolution has been achieved for 
the whole sky. The overall spatial distribution shows 
a clear wavelength dependence. The intensity distri- 
butions of the emission in shorter wavelengths (65 /im 
& 90 /im) show the concentrated distribution on the 
Galactic plane, indicating a tight correlation with star- 
formation activity by tracing high ISRF regions at and 
around star-formation regions. On the other hand, rel- 
atively intense diffuse emission is found in emission at 
longer wavelengths (140 /im & 160 /im), especially from 
the regions at high Galactic latitudes. Emission from 
lower temperature dust in high latitude cirrus clouds is 
thus well traced by longer wavelength observations and 
the ISM total mass distribution is well traced by these 
observations. 
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Figure 1: All-sky FIR intensity maps observed in the 
AKARI diffuse survey. Four bands are displayed from 
top to bottom: N60 (65 fim), WIDE-S (90 fim), WIDE- 
L (140 fim), and N160 (160 fim). The intensity units 
are in MJy sr _1 . 

Zodiacal emission has been subtracted from the im- 
ages and the residual emission can be noticed in shorter 
wavelength images. The residual emission is attributed 
to the Zodiacal dust bands, as we only subtract a spa- 



tially smooth component of the Zodiacal emission and 
do not subtract the emission with smaller scale struc- 
tures. Taking advantage of the higher spatial resolu- 
tions of the AKARI full-sky images, we will investigate 
these structures down to arc-minute spatial scales to 
model and discriminate them from other astronomical 
emissions (Ootsubo et al. 2012). 

We have primarily calibrated the data using the 
more sensitive slow-scan data (Shirahata et al. 2009), 
as the spatial scan speeds are slow enough so that the 
data have better sensitivity and are virtually not af- 
fected by the transient response of the detector. Ko- 
mugi et al. (2012) estimate the 1-a detection limit of 
the survey data as 12 MJy sr _1 , 2 MJy sr _1 , 7 MJy 
sr" 1 , and 4 MJy sr" 1 for N60, WIDE-S, WIDE-L, and 
N160 bands, respectively. The efficacy of the linearity 
correction of the detector signal was also checked with 
slow scan images that observed the same position in 
the sky and is not affected by the detector transient re- 
sponse due to its longer integration time. We confirm 
that a good linearity correction has been achieved for 
all the bands, from the faintest detector signals up to 
the bright signals around the Galactic center of > 10 
[GJy sr -1 ], with relative accuracy of ~ 20%. It is 
important to note that the good linearity correction 
achieved for the data also confirms the successful cor- 
rection of the spatial distortion of the data caused by 
the transient response of the detector. The linearity 
and absolute intensity have been checked with other 
IR observations by ISOPHOT (Kessler et al. 1996), 
IRAS/IRIS (Miville-Deschenes & Lagache 2005), and 
COBE/DIRBE (Hauser et al. 1998). Takita et al. 
(2012) concluded that the correlation of AKARI survey 
data with other FIR data confirms the reliable calibra- 
tion. 

4. Discussion 

The aim of our survey observation is twofold, 1) to ob- 
tain the global distribution of ISM with higher spatial 
resolutions, and 2) to study the SEDs and to evaluate 
the colour temperature of LGs with multi-band pho- 
tometry that has good wavelength coverage. In the 
following discussion, we investigate the characteristics 
of the survey data from these two viewpoints. 

4.1. Spatial Power Spectra 

One of the key characteristics of the AKARI FIR sur- 
vey is its large dynamic range that enables simultane- 
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ous recovery of both large scale spatial structures and 
more localised small scale variations of the FIR emis- 
sion. This characteristic can be examined by calculat- 
ing the spatial power spectra of the cirrus distribution 
taken from the AKARI diffuse all sky survey, as the cir- 
rus power spectra can be well represented by power-law 
spectra (see e.g., Miville-Deschenes et al. 2010; Martin 
et al. 2010). 

Miville-Deschenes et al. (2010) derived the 
IRAS/IRIS 100 urn and Herschel/SPIRE 250 fim com- 
bined spatial power spectrum at the Polaris flare region 
and confirmed that the power-law nature of the cirrus 
power spectrum is kept down to sub-arcminute scales, 
with a power-law index 7 = —2.65 ± 0.10 on scales 
0.025 < k < 2 [arcmin^ 1 ] (also see Martin et al. 2010). 
Since the Polaris flare is a high Galactic latitude cirrus 
cloud with virtually no sign of star-formation activity 
(Martin et al. 2010), the region is well suited to assess 
the nature of the cirrus power spectrum. 

The AKARI FIR image of the Polaris flare region 
and its spatial power spectra are shown in Figs. [2] 
& El The cirrus spectrum can be represented by 




Figure 2: WIDE-L 140/im intensity image of the Polaris 
Flare region in Galactic coordinates centering at I = 
121° and b = 28°. The size of the images are 15° x 15°. 
The intensity units are in MJy sr . 

a single power-law, as indicated by many other au- 
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Figure 3: Spatial power spectra of AKARI intensities 
observed in the Polaris flare region shown in Fig. [5] 
(solid lines). Dotted lines show point spread functions 
(PSFs) of N60 and N160 bands. Higher spatial fre- 
quency components of N60 & N160 power spectra are 
dominated by high frequency noise and well represented 
by their PSFs. 

thors (Miville-Deschenes et al. 2010 and the references 
therein). The power law indices are identical for the 
four wavebands. It is estimated as 7i4Q/«n — 2.55±0.01 
by fitting the 140 fim spectrum in k = 0.01-0.33 
[arcmin -1 ] wavenumber range, and the estimated 7 is 
consistent with that by Miville-Deschenes et al. (2010) 
within the range of error. The spectra clearly show the 
wide dynamic range of our observation, as the spatial 
power component is well retrieved from > 10° large 
scale distribution to ~ 3' small scale structures at 140 
/im and even smaller structures at 90 /im. 

On the other hand, excess components in the power 
spectrum are found at the high spatial frequencies in 
the two narrow bands (N60 & N160). This is due to 
the higher noise levels of these two bands compared to 
the wide bands (WIDE-S & WIDE-L). The spatial dy- 
namic range of the narrow bands is thus limited for the 
faint cirrus component. Better recovery of these high 
frequency spatial component of fainter regions could be 
obtained by adjusting the correction function of the de- 
tector transient response, as the correction is primarily 
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an amplification of the high spatial frequency compo- 
nent and thus causes excess noise at fainter regions (see 
Doi et al. 2009). This adjustment will be considered in 
a future update of the image data. 

4.2. Spectral Energy Distribution 
Another key aspect of the AKAR1 FIR survey is its res- 
olution of the spectral peak of the dust SED with four 
wavebands. Figs. [4] & [5] show the intensity distribution 
of the Cygnus X region and its SED at the centre of 
the region (I — 80°, b = +1°). Ancillary data at the 
same position by IRAS/IRIS (Miville-Deschenes & La- 
gache 2005) and COBE/DIRBE (Hauser et al. 1998) 
and their spatial variation within the 0.7° DIRBE beam 
are also shown. The AKARI data have the largest spa- 
tial variation because of higher spatial resolutions. 

Good coincidence between AKARI and the ancillary 
data is found. The intensity at 140/um is the highest, 
confirming that the peak of dust SED is well covered 
by the AKARI data. WIDE-S, WIDE-L, N160 and 
ancillary data at > 90^tm are well fit with a modified 
blackbody of T= 20.5 [K] and /3 = -2, indicating that 
the AKARI TIR data can be a good tracer of the tem- 
perature and the total amount of LGs. Since AKARI 
has higher spatial resolution comparing to other all- 
sky survey observations as well as huge spatial dynamic 
range at this wavelength range, we can make a detailed 
investigation of the spatial distribution of temperature, 
and the total amount of LGs. 

On the other hand, the N60 observations, and other 
data at 60/im show significant excess from a modified 
blackbody, indicating a significant contribution from 
either higher temperature LGs along the line of sight 
or stochastically heated small grains (SGs). A fit using 
the DustEM SED model by Compiegne et al. (2011) 
with an assumption of uniform incident radiation is also 
shown in Fig. [5] A significant contribution to the ratio 
of emission from SGs to the total observed intensity is 
estimated from the model, not only for N60 (38%), but 
also for WIDE-S (20%). Hence it is important to note 
that the 90//m/140/im colour temperature could over- 
estimate the dust temperature, though the ambiguity 
is much smaller than the widely referred to issues with 
the 60/im/100/im IRAS colour temperature. 

The shorter wavelength excess emission is, how- 
ever, not completely explained by the modelled emis- 
sion from small dust particles with uniform irradiation. 
Temperature variation of dust particles, or intensity 




Figure 4: An intensity map the Cyngus X region in 
Galactic coordinates centering at I = 80° and b = +1° 
with the image size of 15° x 15°. The image is a three- 
colour composite of N60 (blue), WIDE-S (green), and 
WIDE-L (red). 

variations of interstellar radiation field, along the line of 
sight also should be taken into account. The wealth of 
spatial structures and colour variations shown in Fig. U 
clearly indicate the complexity of the ISM and its irra- 
diation environment. The newly achieved AKARI FIR 
high-spatial resolution images of the whole sky should 
be a powerful tool to investigate the detailed spatial 
structure of ISM and its physical environment. 

5. Conclusion 

We have processed a full sky image of the FIR AKARI 
survey at 65/im, 90/im, 140/im and 160/im. The im- 
ages achieve a detection limit of < 10 MJy sr™ 1 with 
absolute and relative photometric accuracies of < 20%. 
Together with the > 99% coverage of the whole sky, 
the high spatial resolution of the FIR AKARI survey 
from 1 to 2 arcmin achieves an unprecedented spatial 
dynamic range of > 10 4 . A good wavelength cover- 
age from 50fim to 180/im with four photometric bands 
makes the data a good tracer of temperature and the 
total amount of the dust component in the ISM. The 
wealth of spatial structures and SED variations found 
in the image clearly indicate the complexity of the ISM 
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Figure 5: Spectrum Energy Distribution at Cygnus 
X region (/ = 80°, b = +1°) in units of AnvI y N n 
ergs s^ 1 H-atonT 1 . The AKARI FIR data as well as 
IRAS/IRIS (Miville-Deschenes & Lagache 2005) and 
COBE/DIRBE (Hauser et al. 1998) data are shown. 
The data are averaged over the 0.7° DIRBE beam, and 
mean values and those standard deviations are indi- 
cated. The red dotted line is an estimation of the LGs 
colour temperature by fitting 90^m, 140/im, 160/im 
data with a single modified blackbody spectrum. The 
blue solid line is a model fit of the four ^4KAis!/data by 
DustEM SED model (Compiegne et al. 2011). The de- 
composition of the model emission into three dust com- 
ponents (PAHs, SGs and LGs) are indicated as black 
dotted lines. Relative contribution of emissions from 
SGs to each AKARI intensity is also indicated in the 
figure. 

and its irradiation environment. The FIR AKARI im- 
age is a new powerful tool to investigate the detailed 
nature of ISM from small scales to the full sky. 
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